NASA Contractor Report 178240 
ICASE REPORT NO. 87 8 

ICASE 


RKH SPACE APPROXIMATIONS FOR THE FEEDBACK 
OPERATOR IN A LINEAR HEREDITARY CONTROL SYSTEM 


J. A. Reneke 
R. E. Fennell 


Contract No. NAS1-18107 
February 1987 


INSTITUTE FOR COMPUTER APPLICATIONS IN SCIENCE AND ENGINEERING 
NASA Langley Research Center, Hampton, Virginia 23665 

Operated by the Universities Space Research Association 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 



RKH Space Approximations for the Feedback 
Operator in a Linear Hereditary Control System* 

J. A. Reneke and R. E. Fennellf 
Department of Mathematical Sciences 
Clemson University 
Clemson S. C. 29634-1907 


Abstract 

Computational implementation of feedback control laws for linear hereditary systems requires the 
approximation of infinite dimensional feedback operators with finite dimensional operators. The 
dense subspaces of K-polygonal functions in reproducing kernel Hilbert spaces, RKH spaces, 
suggest finite dimensional approximations of the matrix representations of the control operators. A 
convergence theorem is developed for the approximations and the numerical implementation of the 
approximations is discussed. 


* This research was supported in part by the following grants AFOSR-84-02336 and 
ONR-N00014-86-K-0693. 

t Partially supported by NASA Contract No. NAS 1-18 107 while the second author was in 
residence at ICASE, NASA Langley Research Center, Hampton, VA 23665-5225. 


i 



RKH Space Approximations for the Feedback 
Operator in a Linear Hereditary Control System 


1. Introduction . Linear hereditary systems are complicated by the fact that state space formulations 
require an infinite dimensional state space. For linear functional differential equations the usual 
choice of state space leads to a description by an abstract evolution equation and to approximations 
by semigroup techniques [2, 5, 8, 9]. We have advocated in several papers an alternate approach, 
i.e., something other than ordinary differential equations in infinite dimensional spaces [15, 16]. 
Our approach is to describe system dynamics by operator equations defined on a reproducing 
kernel Hilbert (RKH) space. This approach, motivated by an integral equation or an input-output 
operator description of die system dynamics is general enough to include linear functional 
differential equations and Volterra integral equations. In this approach one gives up the familiar 
guide provided by finite dimensional state space systems; but one gains the powerful methods 
associated with RKH spaces [1, 10, 13, 18]. In particular, matrix representations of continuous 
linear transformations and the approximations suggested by the convenient dense subspaces of 
K-polygonal functions are available in RKH spaces. 

In this paper we apply these methods to the problem of constructing suitable finite dimensional 
approximations to the feedback operator of a linear hereditary system. We are concerned with 
systems as diagrammed in Figure 1, where A denotes the system input-output operator, D denotes 
a feedback control operator, and f and h denote, respectively, the system input and response. One 
should think of the "physical plant" as being modelled by the input-output operator A. We assume 
A is linear and causal but is not necessarily time invariant or memoryless. 

In applications the feedback operator is determined by a control system designer to achieve 
some desired objective, i.e., stabilizing the system, tracking a desired trajectory, or obtaining a 
desired terminal condition for the system response. Such feedback operators may simply involve 
proportional/integral/derivative terms or more complex dynamic compensation. Synthesis methods 
for hereditary systems [3, 4, 8, 16] lead to feedback operators D which are linear and causal. 
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Figure 1. Feedback Control System. 

In this paper we are not concerned with the synthesis of a feeback operator to achieve a desired 
objective. Motivated by problems of digital implementation of such control operators, we present 
methods to obtain finite dimensional (finite dimensional range) operator approximations to the 
system feedback operator. Digital implementation requires storage of a discrete representation of 
such approximations and rules for their evaluation. In Section 3 projection methods are used to 
obtain feedback operator approximations and show convergence of these and related system 
operator approximations as the sampling interval goes to zero. 

In order that these approximations satisfy the physical requirements of causality a simple 
"hold" of one sampling period is introduced in their definition. Such delays in digital control 
systems result from computational and data processing considerations which can become a 
significant factor in the control of high order systems [12]. Although one normally expects to use 
equally spaced samples, this is not a requirement for our methods. Convergence of related operator 
approximations has been considered in our previous effort [15]. 

In Section 4 we discuss the computational implementation of our approximation scheme. A 
matrix representation of the approximate feedback operator is presented. If the feedback operator is 
time invariant, it is shown that evaluation of the approximate feedback operator can be achieved 
through a simple convolution formula and computer storage requirements can be greatly reduced. 
We illustrate the effect of implementing our feedback operator approximations by considering two 
numerical examples in Section 5. These examples indicate a slight deterioration in the control 
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system performance as the holding time increases and point to several interesting problems for 
future research. 

For finite dimensional state space systems, use of numerical integration rules and related 
approximation schemes to obtain digital filters is well understood [7]. The significance of our 
results lies in the fact that we have provided an analysis of a much broader class of feedback 
operator approximations and a compatible scheme for their digital representation and 
implementation. Examples show that our methods reduce to standard approximations in the state 
space case. 


2. Preliminaries . We begin by presenting our abstract setting. One should refer to the text [19] by 
Willems for a basic discussion of the operator equation approach. For completeness, we present 
only necessary background information. Let r be a nonnegative number, S the number interval [-r, 
«>), d a positive integer, X = R d , <•, •> the usual inner product on X, and | • | the corresponding 
inner product norm. Let G be the class of functions from S into X to which f belongs provided 1) f 
is continuous on [0, °°) and 2) if r > 0 then f is continuous on [-r, 0) and f(0-) exists. Let {N x , x 
in S} be the family of pseudonorms defined on G by N x (f) = lub{ | f(t) |, -r < t < x}, f in G and x 
in S. 


Let B denote the class of linear transformations of G to which B belongs only in case 


1) [Bf](x) = 0 for each f in G and x in [-r, 0] and 

2) for each compact subinterval [u, v] of [0, °°) there is a number b such that 

| [Bf](t) - [Bf](s) | < b / s N x (f)dx, for each f in G and subinterval [s, t] of [u, v]. 

The nondecreasing function bl is called a variation function for B on [u, v]. 

We will reserve C to denote the element of B defined by 


[Cf](t) = 



-r < t < 0 
0<t 
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for f in G and t in S. The simple state space system given formally by h(t) = f(t) + a/ n h(x)dx, 


where a is a d x d matrix and f is in G, would be written as h = f + aCh. The general class of 
hereditary systems to be considered are of the form h = f + Bh with B in B. One should realize that 
our abstract setting is based upon an integral equation description of the system dynamics. 

We can obtain an equivalent class of models as follows. Let A denote the class of linear 
transformations of G to which A belongs only in case A - 1 is in B where I is the identity on G. If 
B is in B then I - B is a reversible map from G onto G and (I - B)' 1 is in A. Similarly, an element 
A of A is a reversible map from G onto G and I - A' 1 is in B. Thus we have equivalent system 
descriptions in terms of B = I - A' 1 or A = (I - B)' 1 [13]. For the system as diagrammed in Figure 
1, if B = I - A' 1 andD is in B then h = A(f + Dh) = (I - B) _1 (f + Dh). Thus h = f + Bh + Dh or h = 
(I - B - D) _1 f. Here (I - B - D)' 1 exists since B + D belongs to B. 

Our analysis takes place in a reproducing kernel Hilbert space associated with the operators in 
A and B. Let k denote the increasing function defined on S by 

1 + r + t if -r < t < 0 

1 + 2r + t if 0 < t 


k(t) = 


{ 


for t in S. Let G H denote the subspace of functions in G which are Hellinger integrable with 
respect to k, i.e., f is in G H only in case for each compact subinterval [a, b] of S there is a number 
M such that 

n 

X i f(s p ) - f<Vi> I 2 1 « s p> - k <Vi” s X 1 df <Vi’ V |2 ' V s M 

p=l s 

for each partition {s p }° of [a, b]. Here df(s pl , s p ) = f(s p ) - f(s p 1 ). The least such number M is 
denoted by jj df | 2 / dk. 

Note that elements of G H are absolutely continuous on compact subintervals of [0, H and 
[-r, 0) if r > 0 [10]. Furthermore, if each of f and g is in G H and [a, b] is a compact subinterval of S 

then the integral f <df, dg> / dk exists as a limit through refinement of partitions s of [a, b] of 
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approximating sums of the form X s <df(s p l , s p ), dg(s p . 1 , s p )> / dk(s p l , s p ). Let G M denote the 

subspace of G H to which f belongs only in case J r | df | / dk < °° and the inner product for 
given by Qjtf, g) = <f(-r), g(-r)> + J r <df, dg> / dk. The space (G^, Q^} is a complete 

inner product space with norm and has a reproducing kernel K given by K(s, t)x = k(min(s, 
t))x for (s, t) in S x S and x in X. Thus K( , t)x belongs to G^ for each t in S and x in X and 
<f(t), x> = Q H (f, K( ,t)x) for each f in G^, t in S and x in X. 

It may be shown that the norm restricted to P T G H T > 0 is equivalent to the Sobolev norm 
f T 2 2 1/2 

(J 0 { I f(0 I +1 f'(t) | } dt) . However use of the norm N m or equivalently the kernel K results 

in simplified representations and approximations of continuous linear operators. 

One can show [15] that an element B of B maps G into G H and hence an element A of A 
maps G h onto G H . We will be concerned with the restrictions of elements of A and B to G H but 
will not introduce any special notation for these restrictions. 

Let {P x , x in S} denote the family of projections on G given by 

{ f(t) -r < t < x 

f(x) x < t 

for each f in G and (x, t) in S x S. It is important to note that operators B in B are causal, i.e., 
P t BP t for each T > 0. 

To illustrate the properties of the kernel K let us compute the adjoint in (G^, Q^} of P T C, 

T > 0. If f is in G M , t is in S, and x is in X then 

<[(P T C)*f](t), x> = Qeo((P T C)*f> K( , t)x) 

= Qoo(f,P T CK (»t) x ) 
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= f <df, dCK( ,t)x> / dk 

* o 

= f <df, K( ,t)x> 

J o 


= <f(T), K(T, t)x> 


- <f(0), K(0, t)x> 


Jo 


dK( ,t)x> 


If -r < t < 0 then <[P T C)*f](t),x> = <f(T) - f(0), k(t)x>. If 0 < t < T then <[(P x C)*f](t),x> = 
<f(T), k(t)x> - <f(0), k(0)x> - <[Cf](t), x>. If T < t then <[(P T C)*f](t), x> = <f(T), k(T)x> - 
<f(0), k(0)x> - <[Cf](T), x>. 

Thus 

[(P T C)*f](t) = K(t, T)f(T) - K(t, 0)f(0) - [P T Cf](t) . 

This expression will be used in the following section. 


3. Matrix representations and the convergence theorem . One of the many advantages of an RKH 
space setting for our analysis is the tractability of the dense linear subspace defined in terms of the 
kernel K, the K-polygonal functions [10]. Consider the family of projections {IT t } from G into 
with t = {t p }° and increasing sequence in S defined by 

{ f(t Q ) -r < s < t Q 

(i/dk(t p _i, t p )){dk( S , yfa^) + dk(t pl , s)f(t p )> t pl < s < t p 
t(t n ) t n < s 

for each f in G and s in S. Elements of II t G are called K-polygons. One can show that for each 
T > -r, the union of the subspaces ri t G H , {t p }° a partition of [-r, Tj is dense in P T G H with 

respect to N^. 

Suppose that D is an element of B, {t p }” is a partition of some interval [-r, T], and 
c = mesh(t) = max (t p - t p l ), p= 1,2, ..., n. When T > 0 we assume that 0 is the partition. Let D t 
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-r < u < c - r 


denote the linear transformation of G defined by 


[D t f](u) = 



0 

[n t Dn t f](u- C ) 


c-r<u 


for each f in G. At first it is natural to choose n^Ilt to approximate D [15] but ntDri t does not 
satisfy the necessary causality requirement. Here D t is causal and we show D t approximates D. 
Thinking of D as a feedback operator and II t f as sampled data, Djf represents an approximation to 
Df obtained from a translation of 11^11^. 

For linear functional differential equations the basic feedback operator arising from linear 

quadratic regulator design [3, 4, 6, 8 ] has the form W(t) = -K^tjhft) - / tr K 2 (t, x)h(x)dx where 

Kj and K 2 denote matrix valued functions. Since our approach is based upon an integral equation 
description of the system dynamics the feedback operator D would be of the form [Dh](t) = 

ft fX 

- J Q (Kj(x)h(x) + r K 2 (u, x)h(u)du} dx. With appropriate boundedness assumptions on K 1 
and K 2 , D belongs to B. 

To further illustrate the possibilities for the feedback operator D, we consider a finite 
dimensional state space example. Robust control system design methods for a plant x p ' = a p x p + 
B p u p , z = y p x p lead to a compensator x c ' = a c x c + B c z, y = y c x c + 8 c z and interconnection 
u p = r rcf - y where r ref denotes a reference signal. Integration yields x p = x p (0) + a p Cx p + B p Cu p 
or x p = (I - a p C) _1 (x p (0) + B p Cu p ), x c = (I - a c CJ’^Cz, and y = {y c (I - a c CY 1 B c C + 5 c }y p x p . 
Letting u = B p Cu p , we obtain u = B p Cr ref - DXp with D = B p C{y c (I - a c C) _1 B c C + 8 C } y p . Here D 
belongs to B and our approximation metohds apply. 


Theorem 1 . The operator D t is in B . Furthermore, any variation function bk of D is also a 
variation function of D t . 
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Proof. If bk is a variation function of D, f is in G, and 0 < Vi* u-c<v-c<t p then 

PVK v) - [D t f](u) | = (l/dk(t p l , t p )) | dk(tp_p v - 0 )^^^) + dk(v - c, t p )[Dn t f](t p . 1 ) - 

dkO^, u ' c)[Dn t f](t p ) - dk(u - c, t p )[Dn t f](t p . 1 ) | 

N s (n,f)ds 

dkf Vi V Vu 

< b dk(u, v) N t(p) (n t f) 

<b dk(u, v) N t(p) (f) 

<b f N(f)dk(s) . 


If 0 < t p .j ^ u - c < tp < tq.j < v - c < t q then 


[DjfKv) - PflCu) | 

^ I [Djf]( V ) - [DjfKtq.! + C) 




+ C) - [D t f](t. x + c) 


+ |[D t f](t p + c)-[D t f](u)| 

<b dk(u, t p + c)N t(p) (f) 

+ b dk(tj _j + c, t + c)N t(i) (f) 

i=p+l 

+ b dk(t j + c, vJN,/ df) 


<b N (f)dk(s) 


and we are through. 


We are concerned in this paper with the problem of substituting the finite dimensional 
approximation D t , i.e., D t has finite dimensional range, for a given feedback operator D, see Figure 


8 


1. We begin by studing the convergence of nets D t , t a partition of [-r, T]. In an application the 
feedback operator is chosen to achieve some desired objective, i.e., the plant h = f + Bh + Dh, or h 
= (I - B - D) _1 f with B = I - A' 1 behaves satisfactorily. Since D t represents a digital implementation 
of D, it is important to also consider convergence of the nets (I - B - Dj)' 1 to the system 
imput/output operator (I - B - D)' 1 . In later sections we take up a computational implementation of 
these ideas and apply then to two examples. 

A function L from S x S into the continuous linear transformations of X is called the matrix 
representation of D provided <[Df](t), x> = Q M (f, L( , t)x) for each f in G^, t in S, and x in X. A 
fundamental property of RKH spaces is that every continuous linear transformation between such 
spaces has a matrix representation. A problem is that D maps g h into G h rather than G m into G TO . 

For T in S, P T D is continuous from G ro into G TO with respect to and so P T D has a matrix 
representation. Note that <[P T Df](t), x> = Q M (P T DP T f, K( , t)x) = Q oo (P T f, (P T D)*K( , t)x), for 
-r < t < T, f in G H , and x in X. Let L(s, t) x = [(P T D)*K( , t)x](s) for -r < s, t < T and x in X. To 
obtain L without having to compute (P T D)* note that, for 1 < i, j < d and -r < s, t < T, < e ; , L(s, 
t)ej > = Q oc (P t DK( , s)e i; K( , t)e-) = <[DK( , sje^t), e.>, where e ; and e^ are standard basis 
elements of X = R d . In light of this equation we speak of the matrix representation of D rather than 
of P t D. 

Note that L(s, t) = 0 for -r < t < 0. For L causality means that given -r < t < s, <x, L(s, t)y> = 
<[DK( ,s)x](t),y> = <[DK( ,t)x](t),y> = <x, L(t, t)y> for each x and y in X, i.e., L(s, t) = L(t, t). 
Hence L( ,t)x is in G^ for each t in S and x in X. 

For each positive number T, let N^ denote the operator pseudonorm defined on the linear 
operators F of G H which are continuous with respect to N T by N HT (F) = sup{N T (F g )/N oo (P Tg ) | g in 
G H .P Tg * 0}. Assume for the rest of the paper that T is a fixed positive number. We will discuss 
convergence of D t and the control systems {(I - B - D t ) -1 }, t a partition of [-r, T], using N^ and 
the matrix representations L and M of D and D t , respectively. 
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Theorem 2 . The operator P T D is the limit through refinement of the net of finite dimensional 
operators {D t }, t a partition of [-r, T], with respect to the operator pseudonorm N HT . 

Lemma . If t is a partition of [-r, T], u is in the range of t, and x is in X then N m (L( ,u)x - 
M( , u + c)x) 2 = ,u)x) 2 - X t | dL( , u) x | 2 /dk. 

Proof . If y is in X, v is in S, and c = mesh(t) then 

<y, M(v, u + c)x> = <[D t K( , v)y](u + c), x> 

= Q J]1 t Dn t K( , v )y, K(,u)x) 

= QJ K ( > v)y, IlfLC , u)x) 

= <y, [IIjLC , u)x](v)> . 

Thus 

Q m (L( , u)x, M( , u + c)x) = QooCLC , u)x, n,L( , u)x) 

= Z t I dL( , u)x [ 2 / dk 
= N to (M( , u + c)x) 2 . 

Hence N^CLC , u)x - M( , u + c)x) 2 = NJLC , u)x) 2 - Z t | dL( , u)x | 2 / dk . 

L emm a. If [u, v] is a subinterval of [0, T] and x is in X then N m (L( , v)x - L( , u)x) < 
N^CPjD) | x | dk(u, v) 1/2 . Also, for each partition t of [-r, T], N oo (M( , v)x - M( , u)x) < 
N TO (P T D)|x|dk(u,v) 1/2 . 

Proof. Note that NJL( , v)x-L( , u)x) = NJ(P T D)*(K( , v) - K( , u))x) < NJPjD) | x |- 
dk(u,v) 1/2 and , v)x - M( , u)x) = N oo (IT t (P T D)*n t (K( , v - c) - K( , u - c))x) < 

N m (P t D) | x | dk(u, v) 1/2 . 
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Proof of the theorem. If {s p }™ is a partition of [-r, T], {t q }° is a refinement of s, c = mesh(t), 
f is in G h , x is in X, and u is in [s pl , s p ] for some positive integer p then 

| <[Df](u) - [D/Ku), x> | 

= I Qoo(f> L( , u)x - M( , u)x) I 
^ N oo(P T 0 (NJL( , u)x - L( , s p . 1 )x) 

+ N to (L( , Sp.^x - M( , Sp.! + c)x) + N..(M( , s p l + c)x - M( , u)x)} 
<NJP T f){NJP T D) | x | dk(s p4 , u) 1/2 + (NJL( , s^x) 2 

- I t | dL( , s^j) x | 2 /dk) 1/2 + NJP T D) | x | (dk^, u) + dk(0, c)) 1/2 } 

* NJP T f){3N 00 (P T D) | x | (mesh(s)) 1/2 + (N^(L( , s^x) 2 

- Z t | dL ( , s p _j) x | 2 /dk) 1/2 } . 

Suppose that e > 0, {s p }g is a partition of [-r, T] such that mesh(s) < (£76(1^ (P^D) + l)) 2 /d, 

t' is a refinement of s such that if t refines t', p = 1, 2, ..., n and i = 1,2, d then 
N 00 (L( , s p . 1 ))e i ) 2 - £ t | dL( , Sp.^ej | 2 / dk < (e/2) 2 /d. If t refines t' and | x | = 1 then 
| < [Df](u) - [D t f](u), x> | 

d 

^ X I x i 1 1 < t Df K u ) “ PtHCu), e t > | 

i=l 
d 

s ^ I Xi | N„CP T 0 e /d 1/2 

i=l 

^elxIN^f) , 
i.e., N h T (D - D t ) < e. 

Theorem 3 . If B is in B then (I - B - D)' 1 is the limit through refinement of the net {(I 
t a partition of [-r, T], with respect to the operator pseudonorm N^. 

Lemma . If bk is a variation function for D then be bT k is a variation function for (I 

KT]. 


-D-D^ 1 }, 
- D)’ 1 - 1 on 
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Proof . If h = f + Dh then by the Gronwall inequality N u (h) < e bu N u (f) for each u in [0, T], 
Thus for each subinterval [u, v] of [0, T] 


[(I - D)-!f](v) - f(v) - [(I - D) _1 f](u) + f(u) | 
<b N ((I - D^dkCs) 

•'ll 

<b e bs N(f)dk(s) 

•'ll 

f V 

<be bT N s (f)dk(s) . 


Proof of Ihg theorem . Assume bjk and b 2 k are variation functions for B and D, respectively. 
Then N T ((I - B - D t ) -1 ) < 1 + be bT (k(T) - k(0)), where b = bj + b 2 . Hence for each f in G H . 

N t ((I - B - D t ) _1 f - (I - B - D) _1 f) 

< N t ((I - B - - D t )N M (P T (I - B - D)' 1 )N 00 (P T f) 

< {1 + be bT (k(T) - k(0))}N HT (D - D t )NJP T (I - B - D)' 1 ) N M (P T f), 

•i.e., NhxCCI - B - D t )- J - (I - B - D)' 1 ) < {1 + be bT (k(T) - kCOWN^I - B - D)‘ 1 )N HT (D - D t ). 
Thus (I - B - D)' 1 is the limit with respect to the operator pseudonorm N m of the net {(I - B - 
Dj)' 1 }, t a partition of [0, T]. 

Thus N t (I - B - D t ) _1 f - (I - B - D) _1 f) converges to zero with respect to refinements for each f 
in Gjjj i.e., the solution of h = f + Bh + D t h converges uniformly on [-r, T] to the solution of h = f 
+ Bh + Dh. 

4. Computational implementation . In this section we elaborate on the discrete structure of the 
approximation D t which facilitates a digital implementation. Recall that L and M denote the matrix 
representations of D and D t respectively. In this paper we are not concerned with computing L(s, t) 
or M(s, t) for various values of s and t in S but with the effect of using the approximations D t in 
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place of D. We would expect that most D’s met in practice would have representations in terms of 
familiar integrals and we can obtain the various values of L and M using some existing numerical 
quadrature code. 

In addition, the extensive literature [2, 3, 8, 9] on numerical approximations for hereditary 
systems is available for evaluations of those L and M arising from linear quadratic regulator 
designs. 

If (t }” is an increasing sequence in S, t . < s < t , and f is in G then [IlD]Tf](s) = 

(dk(s, t p )[Dll t f](tp. 1 ) + dk(t pl , s)[Dn t f] (t p ))/dk(t p _ 1 , t p ). Thus to evaluate 11,01^ we need to 
know [Dll t f[(tp) for p = 0, 1, 2, ..., n. Let f t denote the (n + l)d vector given by f t = col(f(t 0 ), 
f(tj), ..., f(t n )) and L t the (n + l)d x (n + l)d block matrix whose (p, q) block is given by L t (p, q)^ 

= <L(t p , t q )ej, e L > for 0 < p, q < n and 1 < i, j < d. Define K t in an analogous way. 

We have IT t f = L q=0 K( , t q )c q where c q = [(K t ) _1 f t ](q) for q = 0, 1, 2, ..., n. Hence 

Dn t f = X q=0 DK( , t q ) Cq = X q=0 L(t q , )*c q and, for p = 0, 1, ..., n [D]l t f](t p ) = [(L t )*(K t )\](p). 

Let M t denote the (n + l)d x (n +l)d block matrix defined by 

{ 0 P = 0 

-l 

[(L t )*(K t ) j(p - 1, q) 1 <p < n, 0<q<n 

Then [D t f](t p ) = [M t f t ](p). 

The matrices (L t )*(K t ) _1 and M t denote discrete representations of the operators Il t Dn t and D t . 
respectively. Let us look closer at the structure of these matrices. 

One can show that (K t ) _1 is symmetric block tridiagonal. For p = 1,2, ..., n. the p-th block in 
the first subdiagonal is -(l/dk(t p l , t ))I. The first main diagonal block is (l/k(t 0 ) + l/dk(t 0 , tj))I, 
the last main diagonal block is (l/dk(t n _p t n ))I, and, for p = 2, ..., n the p-th main diagonal block is 
(l/dk(t p . 2 , tp.^ + l/dk(t pl , t p ))I. If 0 <p < q<n then [(L^K^Kp, q) = (L t (p, p))*{[(K t )- 1 ](q -1, 
q) + [(K^Kq, q) + [(K^Kq + 1), q)> = 0, i.e., (L t )*(K t )- ! is block lower triangular. 

Clearly M t is also block lower triangular. 
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If D is a time invariant operator, i.e., [Dg](t + b) = [Df](t) with g(t) = f(t - b) for each f in (I - 
P 0 )G h , b > 0, t > 0 then the structure of the matrices simplifies further. In this case L(v, w) - L(u, 
w) = L(v + b, w + b) - L(a + b), w + b) for 0 < u < v < w and b > 0 [15]. Let t p = - r + pc, p = 0, 


..., n with r = p'c. It follows that 
(L t )*(K t ) 1 (p,q) = 0 

= {L t (p'+1, p +1)* - L t (p', p'+l)}/(r+c) 

= {L t (p +1, p'+l)* - L t (p, p +l)*}/c 
= {L t (p'+1, p)* - L t (p y , p)*}/(r+c) 

’ (L t (p +1> P‘l)* ■ L t (p', p-l)*}/c 
= {L t (p +1, p'+l+j)* - L t (p', p'+l+j)*}/c 
-{L t (p'+l,p , +j)*-L t (p',p , +j)*}/c . 

Thus the storage requirements for (L t )*K t ) _1 and M t are significantly reduced. One need only 
store (L t )*(K t ) _1 (p,q) for p = p' < q < n, p = p' + 1 < q < n, 0 < q < p' < p. Here [DjfKtp) = 0 for 


0 < q < p < p' 
p = q = p' + 1 

p = q > p + 1 

p > q = p + i 
p = q+j>q>p 


p < p' + 1 and for p' + 1 < p < n 


[D,f](t p ) = (L,)*(K,)' l (p-l, q)f(t q ) 

1 ^ip-l-(p'+2) 

* 2 q . 0 c-m' <p- >• + a.„ a /w 


where aj = ((L t )*(K t ) _1 (n, n-j) for j = 0, ..., n - (p'+l). This expression simplifies further if r = 0, 
which would be the case if D is of the form [Df](t) = J 0 g(t - s)f(s) ds for example. If r = 0 then 


[ D t f KV = 2, j=0 a j f((p - 1 - j)c) 

for p' + 1 < p < n. Thus in the time invariant case evaluation of [D t f](t p ) involves a convolution 
sum and, consequently, computation time can be reduced through use of an appropriate fast 
algorithm. 
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In the examples to follow, we solve the system h = ]~I t 'f + (B + D t ) t - h with t' a refinement of t 
to obtain approximate solutions of the system h = f + Bh + D t h. 

5. Numerical examples . In order to simulate h = f + Bh + D t h with t = {t p }° a partition of [-r, T] 

one should choose, depending on the nature of B, a differential equation solver, see for example 
[12]. On the interval [-r, 0] we have h = f. Let y = 0 < y +1 . On the interval [0, y +1 ] we know 
f and D t h and so can solve for h. Proceeding in this manner one steps across the interval [0, Tj. 

The following examples, were chosen to illustrate the effect of implementing the 
approximations D t in a feedback control system. Because the examples are simple we chose to 
base our calculations on the corollary to Theorem 3. 

Example 1. Consider the time invariant system h = f + Bh + Dh with [Bh](t) = -4 ^ (t-x)h(t) dx + 

J 0 h(x-.25) dt and [Dh(t) = -4 J 0 h(x) dx for t > 0. Here the feedback operator D stabilizes the 

system h = f + Bh. The effect of this stabilizing feedback is illustrated in Figure 2. In Figure 2 h 
denotes the approximate solution of the system with D = 0 and h 2 , h g , h 32 denote approximate 
solutions of the feedback system with c = t p - yj = 1/128, 1/32, 1/8 respectively. It is interesting 
to note that stability is maintained for relatively large values of the "holding time" c. Here we 
obtain N 2 (fi-h p ) < .031036, .061612, .12045, .22959, .41532 for c = p/256, p = 2, 4, 8, 16, 32 

respectively, where fi denotes the solution of h = f + Bh + Dh. 


Example 2. Let [Bh](t) = -4|| ) (t-x)h(x) dx - J^h(x-.25) dx. For a given disturbance f, we wish to 

choose a control u to drive the response of the system h = f + Bh + u to zero over a terminal 
interval [1.25, 1.5]. Optimization techniques [14] can be used to obtain such a control in the form 
u = Dh + g, i.e., u is the sum of a feedback term and an open loop term. Here D = A _1 C 2 with A 1 
= I - B and g depends on the known disturbance f and known parameters. The effects of 
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More specific knowledge of the feedback operator leads to further simplifications. For 
example, if [Dh](t) = / 0 Kh(x)dx, which is the case for a state space system h'(t) = ah(t) + Bv(t) 

with feedback v(t) = k h(t), then [D t h](c) = 0 and for p > 1 [D t h](pc) = k (c/2){h(0) + 2h(c) + ... + 
2h((p-2)c) + h((p-l)c), i.e., the approximation scheme reduces to the trapezoid rule. 
Implementation of this approximation yields h'(t) = a h(t) + Bk (h((p-l)c) + h(pc)}/2 for pc < t < 
(p+l)c, i.e., a simple "first order hold" approximation for the feedback operator. 

Before illustrating the effect of implementing the approximation D t in a feedback control 
system, we present a corollary to Theorem 3 which provides a method to obtain approximate 
solutions to the system h = f + Dh with D in B. 

For t = {t p }g a partition of [-r, T], if h = ll t f + D f h then h is in Il t G H and for 0 < p < n 

h(t p ) = f(t p ) + [D t h](t p ) or h t (p) = f t (p) + [M t h t ](p). Thus h t = f t + M t h t or h t = ((I - D^IIA = 
^(n+l)d ‘ % 


Corollary to Theorem 3 . If f is in G H , h = (I - D)' 1 f, and e is a positive number then there is a 
partition s of [0, T] such that if {t p }g refines s and 0 < p < n then | h(t p ) - [(I( n+1 ) d - M t ) V t ](p) | < e. 
Proof . Using Theorem 3 with B = 0, there is a partition s of [0, T] such that if {^t}^ refines 

s then Nnj^I - D)' 1 - (I - D t ) _1 ) is smaller than either 1 or (1/2) eKN^if) + 1) and N^f - n t f) < 1/2 
e/CN^CI - D)' 1 ) + 1). Thus 

N-rai-D^f-a-D^ritf) 

< Njn-CCI - D)’ 1 - (I - D^NJf) 

+ N HT ((I-D t )- 1 )NJf-n t f) 

< 1/2 e + N m ((I - D)' 1 + UNJf - 
<e • 

Since ((I - D t ) = (I( n +i)d - M t )f t we are through. 


15 


implementing this control and the approximations u = D t h + g with f(t) = 1 + t for -.25 < t < 1.5 are 
illustrated in Figure 3. Here h denotes the approximate solution of the system with D = 0 and h 2 , 
h 16 , and h 32 denote approximate solutions with c = t p - t pl = 1/128, 1/16. 1/8 respectively. Note 
that the desired terminal condition holds when c = 1/128 and there is a slight deterioration in the 
control system performance for larger values of the "holding time" c. For the example we obtain 
Nj 5 (h - h p ) < .0032262, .0072059, .015469, .032368, .067435, .14337 for c = p/256 and 
p = 2, 4, 8, 16, 32 respectively, where fi denotes the solution of h = f + Bh + Dh + g. 

6. Concluding Remarks . We have presented a theory for the approximation of feedback control 
operators which arise in the analysis of linear hereditary systems. The reproducing kernel Hilbert 
space setting introduced facilitates the analysis and digital representation of such approximations. 
Advances in on board computational capabilities indicate that storage of our discrete representations 
of a feedback operator is feasible. Furthermore, in the design of feedback controllers for infinite 
dimensional systems one need not be restricted to the search for simple feedback gains as in the 
finite dimensional case. 

Although we have not addressed control design issues in this paper, such problems have been 
discussed using our abstract setting within a number of contexts [11, 16]. We are presently 
investigating the robustness of control designs based upon an input/output description of the 
systems. The results and examples presented in this paper indicate that such designs should have 
good robustness properties. 

A unified discussion of our approach to system analysis, approximation, and control will 
appear in the lecture notes, "Structured Hereditary Systems" [17], 
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Figure 2. Approximate solutions for Example I. 
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Figure 3. 
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